How well do we know the electromagnetic form factors of the proton? 
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Several experiments have extracted proton electromagnetic form factors from elastic cross section 
measurements using the Rosenbluth technique. Global analyses of these measurements indicate 
approximate scaling of the electric and magnetic form factors {[i p Ge p /Gm p ~ 1), in contrast to 
recent polarization transfer measurements from Jefferson Lab. We present here a global reanalysis of 
the cross section data aimed at understanding the disagreement between the Rosenbluth extraction 
and the polarization transfer data. We find that the individual cross section measurements are self- 
consistent, and that the new global analysis yields results that are still inconsistent with polarization 
measurements. This discrepancy indicates a fundamental problem in one of the two techniques, or a 
significant error in polarization transfer or cross section measurements. An error in the polarization 
data would imply a large error in the extracted electric form factor, while an error in the cross sections 
implies an uncertainty in the extracted form factors, even if the form factor ratio is measured exactly. 

PACS numbers: PACS number: 25.30.Bf, 13.40.Gp, 14.20.Dh 



I. INTRODUCTION 

The electromagnetic structure of the proton is de- 
scribed by the electric and magnetic form factors. Over 
the past several decades, a large number of experiments 
have measured elastic electron-proton scattering cross 
sections in order to extract the electric and magnetic form 
factors, Ge (Q 2 ) and Gm (Q 2 ) (where — Q 2 is the four- 
momentum transfer squared), using the Rosenbluth tech- 
nique The electric and magnetic form factors have 
been extracted up to Q 2 « 7 GeV 2 by direct Rosenbluth 
separations, and these measurements indicate approxi- 
mate form factor scaling, i.e., (i p Ge /Gm ~ 1 (where 
fi p is the magnetic dipole moment of the proton) , though 
with large uncertainties in Ge p at the highest Q 2 val- 
ues 0,u 

More recently, elastic electron-proton polarization 
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FIG. 1: Ratio of electric to magnetic form factor as extracted 
by Rosenbluth measurements (hollow squares) and from the 
JLab measurements of recoil polarization (solid circles). The 
dashed line is the fit to the polarization transfer data. 



transfer measurements have been performed to obtain 
the ratio Ge p /Gm p - A low Q 2 measurement at MIT- 
Bates obtained values of Ge p /Gm p consistent with 
previous Rosenbluth separations. Later experiments at 
Jefferson Lab (JLab) extended these measurements up to 
Q 2 = 5.6 GeV 2 0,0, 0- and show significant deviations 
from form factor scaling. They show a roughly linear 
decrease of the value of /i p Ge /Gm from unity at low 
Q 2 to approximately 0.3 at Q 2 = 5.6 GeV 2 . Figure fl] 
shows the JLab polarization transfer measurements from 
refs. |U, along with a global Rosenbluth analysis of 
the cross section measurements 0- While the polariza- 
tion transfer technique allows much better measurements 
at high Q 2 values, there is a significant discrepancy even 
in the region where both techniques have comparable un- 
certainties. 

When we combine the cross sections with polarization 
transfer measurements to extract the form factors (see 
Sect. MI Fl or Ref. Q), we find that the values obtained 
for Ge p are significantly different, while Gm p differs only 
at the few percent level, compared to extractions that 
use only the cross sections. Clearly, it is necessary to un- 
derstand the discrepancy in the extracted ratio before we 
can be confident in our knowledge of Ge and, to a lesser 
extent, Gm p - The Rosenbluth data is more sensitive to 
systematic uncertainties, and it has been suggested that 
the different Rosenbluth extractions are inconsistent, and 
thus unreliable. We will examine the consistency of the 
Rosenbluth measurements to test this suggestion. How- 
ever, even if it is demonstrated that the cross sections 
going into the Rosenbluth extractions were incorrect, it 
would not completely solve the problem. The polariza- 
tion measurements determine only the ratio of Ge p to 
Gm p , and so reliable cross sections are still needed to 
extract the actual values of the form factors. Finally, if 
the discrepancy arises from a fundamental problem with 
either of these techniques, it may have implications for 
other measurements. 

The goal of this analysis it to better understand the 
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discrepancy between the Rosenbluth and polarization 
transfer results. We begin by demonstrating that the 
individual Rosenbluth measurements yield consistent re- 
sults when analyzed independently, so that the normal- 
ization uncertainties between different measurements do 
not impact the result. We then perform a global analy- 
sis of the cross section measurements, and determine that 
the results cannot be made to agree with the polarization 
results by excluding an small set of measurements, or by 
making reasonable modifications to the relative normal- 
ization of the various experiments. The paper is orga- 
nized as follows: In section [H] we will review the two 
techniques and summarize the current measurements of 
the form factors. In section Hill we present a new Rosen- 
bluth analysis of the cross section measurements, and 
compare this to the polarization transfer results and ex- 
amine possible scenarios that might explain the discrep- 
ancy between the techniques, such as problems with indi- 
vidual data sets or improper treatment of normalization 
uncertainties when combining cross sections from differ- 
ent experiments. In section HVI we will discuss the results 
of the analysis and implications of the discrepancy be- 
tween the two techniques. Finally, in section we sum- 
marize the results and discuss further tests that can be 
performed to help explain the disagreement between the 
techniques. 



II. OVERVIEW OF FORM FACTOR 
MEASUREMENTS 

We begin with a brief description of the Rosenbluth 
separation and recoil polarization techniques, focusing on 
the existing data and potential problems with the extrac- 
tion techniques. 



A. Rosenbluth Technique 

The unpolarized differential cross section for elastic 
scattering can be written in terms of the cross section 
for scattering from a point charge and the electric and 
magnetic form factors: 
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where 9 e is the electron scattering angle, and r = 
Q 2 /AM 2 . One can then define a reduced cross section, 



da e(l + r) 

d^l CTMott 



= tG 2 (Q 2 ) + eG 2 E (Q 2 ), (2) 



where e is the longitudinal polarization of the virtual pho- 
ton (e- 1 = 1 + 2(1 + t) tan 2 (0 e /2)). At fixed Q 2 , i.e., 
fixed r, the form factors are constant and an depends 
only on e. A Rosenbluth, or longitudinal-transverse (L- 
T), separation involves measuring cross sections at sev- 
eral different beam energies while varying the scattering 



angle to keep Q 2 fixed while varying e. G 2 E can then be 
extracted from the slope of the reduced cross section ver- 
sus e, and tG % m from the intercept. Note that because 
the G 2 M term has a weighting of r/e with respect to the 
G 2 E ^ term, the relative contribution of the electric form 
factor is suppressed at high Q 2 , even for e = 1. 

Because the electric form is extracted from the dif- 
ference of reduced cross section measurements at vari- 
ous e values, the uncertainty in the extracted value of 
G\ (Q 2 ) is roughly the uncertainty in that difference, 

magnified by factors of (Ae) _1 and (tG 2 m /G Ep ). This 
enhancement of the experimental uncertainties can be- 
come quite large when the range of e values covered is 
small or when r (= Q 2 /AM 2 ) is large. This is especially 
important when one combines high-e data from one ex- 
periment with low-e data from another to extract the 
e-dependence of the cross section. In this case, an error 
in the normalization between the data sets will lead to 
an error in G 2 E for all Q 2 values where the data are com- 
bined. If fJ, p GE = Gm , Ge contributes at most 8.3% 
(4.3%) to the total cross section at Q 2 = 5(10) GeV 2 , so a 
normalization difference of 1% between a high-e and low- 
e measurement would change the ratio \i v Ge /Gm by 
12% at Q 2 = 5 GeV 2 and 23% at Q 2 = 10 GeV 2 , more if 
Ae < 1. Therefore, it is vital that one properly account 
for the uncertainty in the relative normalization of the 
data sets when extracting the form factor ratios. The 
decreasing sensitivity to Ge p at large Q 2 values limits 
the range of applicability of Rosenbluth extractions; this 
was the original motivation for the polarization trans- 
fer measurements, whose sensitivity does not decrease as 
rapidly with Q 2 . 



B. Recoil Polarization Technique 

In polarized elastic electron-proton scattering, 
p(e, e'p), the longitudinal (Pi) and transverse (P t ) 
components of the recoil polarization are sensitive to 
different combinations of the electric and magnetic 
elastic form factors. The ratio of the form factors, 
Ge p /Gm f , can be directly related to the components of 



the recoil polarization 
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where E e and E' e are the incoming and scattered elec- 
tron energies, and Pi and Pt are the longitudinal and 
transverse components of the final proton polarization. 
Because Ge /Gm is proportional to the ratio of polar- 
ization components, the measurement does not require 
an accurate knowledge of the beam polarization or ana- 
lyzing power of the recoil polarimetcr. Estimates of ra- 
diative corrections indicate that the effects on the recoil 
polarizations are small and at least partially cancel in the 
ratio of the two polarization component 
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of the polarization transfer data comes from the same 
experimental setup, it is in principle possible that an un- 
accounted for systematic error could cause a false Q 2 - 
dependence in the ratio. There are no known problems 
or inconsistencies in these measurements and this tech- 
nique. At this time, there is no explanation for the differ- 
ent results obtained by the two techniques. If we do not 
understand this discrepancy, then it is difficult to know 
how to correctly combine the polarization transfer mea- 
surements with the cross section measurements in order 
to extract the individual form factors. 



III. REANALYSIS OF THE ROSENBLUTH 
MEASUREMENTS 



FIG. 2: Ratio of electric to magnetic form factor as ex- 
tracted by recoil polarization measurements at MIT-Bates 
(solid squares) and JLab (hollow symbols). The solid line is 
the fit to (diamonds and triangles) from Ref. 8] , while 
the dashed line is a combined fit, including the systematic un- 
certainties (assumed to be fully correlated within each data 
set). 



Figure 13 shows the measured values of u p Ge p /Gm p 
from the MIT-Bates and JLab 00 experi- 
ments, both coincidence and single arm measurements, 
along with the linear fit of Ref. to the data from 00: 

UpGeJGm, = 1 - 0.13(Q 2 - 0.04), (4) 

with Q 2 in GeV 2 . Comparing the data to the fit, the 
total x 2 is is 34.9 for 28 points, including statistical er- 
rors only. Assuming that the systematic uncertainties 
for each experiment are fully correlated, we can vary the 
systematic offset for each data set and the total \ 2 de- 
creases to 33.6. If we allow the systematic offset to vary 
for each data set and refit the Q 2 -dependence to all four 
data sets using the same two-parameter fit as above, i.e., 



/J- P G Ep /G 



M„ 



l-a(Q 2 -Q 2 ), 



(5) 



for Q 2 > Qq, unity for Q 2 < Q\. We obtain a = 
0.135 ± 0.008, Ql = 0.24 ± 0.08, and x 2 =28.1 for 26 de- 
grees of freedom (dashed line in Fig. [SJ. The systematic 
offsets are small (consistent with zero) for all of the data 
sets except the low Q 2 JLab measurement 0, which is 
increased by nearly the full (correlated) systematic un- 
certainty This fit not only has a better x 2 , but also 
decreases the deviation from unity at very low Q 2 val- 
ues, which improves the agreement with the very precise 
Rosenbluth results available below Q 2 = 1.0 GeV 2 . 

The main systematic uncertainties come from inelas- 
tic background processes and determination of the spin- 
precession, both of which have been carefully studied and 
accounted for in the JLab measurements. While this 
technique should be less sensitive to systematic uncer- 
tainties than the Rosenbluth extractions, the discrepancy 
appears at relatively low Q 2 values, where both tech- 
niques give equally precise results. Because almost all 



The global Rosenbluth analysis shown in Fig. ^ may 
disagree with the polarization transfer results for a vari- 
ety of reasons: inclusion of bad data points or data sets in 
the fit, or improper constraints on the relative normaliza- 
tion of data sets. To better understand the discrepancy, 
we have performed a reanalysis of the Rosenbluth mea- 
surements. We will use this to look for errors or inconsis- 
tencies in the data sets and to test possible explanations 
for the discrepancy between the two techniques. 

An initial analysis reproduced the results of the previ- 
ous global fit At this point, several modifications 
were made to the data set for subsequent fits: radiative 
corrections were updated for some of the older measure- 
ments, certain experiments were subdivided into separate 
data sets, some normalization uncertainties were updated 
from those used by Walker , several cross section mea- 
surements were updated with the final published results, 
and a set of data points were excluded. These modifica- 
tions are described in detail in the following section. 



Data Selection 



The global analysis presented here is similar to the one 
presented in Refs. . Table [I] shows the data sets in- 
cluded in the fit, along with a summary of the kinematics 
for each data set. The experiments included are, for the 
most part, the same as in the previous analysis. Two 
additional data sets with measurements in the relevant 
Q 2 region have been included [Trl Il7| . For two of the 
experiments included in the previous fit, we use the final 
published cross sections (T^, , which were not available 
at the time of the previous global analysis. 

For each data set included in the fit, an overall nor- 
malization or scale uncertainty was determined, separate 
from the point-to-point systematic uncertainties. This 
normalization uncertainty is given in, or was estimated 
from, the original publication of the data. In most cases, 
we use the same scale uncertainty as in the previous 
global analysis. For six |U El |H |H H3 of the 
sixteen experiments, the published uncertainties included 
the normalization uncertainties. In the previous fit, these 
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TABLE I: Experiments included in global fit, along with the 
Q 2 range, electron scattering angle range, and normalization 
uncertainty. 



Reference 


Q 2 


e 


Norm. 


Laboratorj 




[GeV 2 l 


[degrees] 


Uncert. 




Janssens-1966 [20] 


0.16-1.17 


40-145 


1.6% 


Mark 111 


Bartel-1966 [21] 


0.39-4.09 


10-25 


2.5% 


DESY 


Albrecht-1966 [22] 


4.08-7.85 


47-76 


8.0% 


DESY 


Albrecht-1967 [23] 


1.95-9.56 


76 


3.0% 


DESY 


Litt-1970 [241 


1.00-3.75 


12-41 


4.0% 


SLAC 


Goitein-1970 1 [25] 


0.27-1.75 


20 


2.0% 


CEA 


Goitein-1970 1 [25] 


2.73-5.84 


19-34 


3.8% 


CEA 


Berger-1971 [261 


0.08-1.95 


25-111 


4.0% 


Bonn 


Price-1971 [27] 


0.27-1.75 


60-90 


1.9% 


CEA 


Bartel-1973 2 [28] 


0.67-3.00 


12-18 


2.1% 


DESY 


Bartel-1973 2 [28] 


0.67-3.00 


86 


2.1% 


DESY 


Bartel-1973 2 [28] 


1.17-3.00 


86-90 


2.1% 


DESY 


Kirk- 1973 [29J 


1.00-9.98 


12-18 


4.0% 


SLAC 


Stein-1975 [16] 


0.10-1.85 


4 


2.4% 


SLAC 


Bosted-1990 [30] 


0.49-1.75 


180 


2.3% 


SLAC 


Rock- 1992 [17] 


2.50-10.0 


10 


4.1% 


SLAC 


Sill- 1993 [18] 


2.88-31.2 


21-33 


3.0% 


SLAC 


Walker-1994 3 [2] 


1.00-3.00 


12-46 


1.9% 


SLAC 


Andivahis-1994 4 [19] 
Andivahis- 1994 s [19] 


1.75-7.00 


13-90 


1.8% 


SLAC 


1.75-8.83 


90 


2.7% 


SLAC 



1 Split into two data sets (see text). 

2 Split into three data sets (see text). 
s Data below 20° excluded. 

4 8 GeV spectrometer. 
° 1.6 GeV spectrometer. 



uncertainties were double counted when an additional 
normalization uncertainty was added. For these experi- 
ments, we apply the same normalization uncertainty, but 
remove it from the published (total) uncertainties to ob- 
tain the point-to-point uncertainties. 

Two of the experiments [l9|, |28| included data taken 
with more than one detector. There will therefore be 
different normalization factor for the data taken in the 
different detectors. In Ref. 0]> these normalization fac- 
tors were measured by taking data at identical kinemat- 
ics for two Q 2 points. We split the experiment into two 
data sets, and fit the normalization factor for each one 
independently. This will allow the normalization factor 
to be determined from both these direct measurements 
and the comparison to the full data set. Because we 
do not apply the normalization factor determined from 
the original analysis, we add a 2% normalization uncer- 
tainty (in quadrature) to the 1.77% uncertainty quoted 
in the original analysis. While this may underestimate 
the uncertainty in the normalization, the result would 
tend to be a larger cross section for this low-e data set, 
which would lead to a smaller value for Ge p /Gm p - As 



will be shown, even with this possible bias towards lower 
values of Ge p /Gm p , the resulting ratio is clearly higher 
then the polarization transfer results. Similarly, the elas- 
tic cross sections in Ref. [28| include three different sets 
of data: electrons detected in a small angle spectrome- 
ter, electrons detected in a large angle spectrometer, and 
protons detected in the small angle spectrometer (corre- 
sponding to large angle electron scattering). We divide 
this experiment into three data sets, each with its own 
normalization factor. Finally, after an initial analysis, it 
was observed that the data from Ref. [2f| was taken un- 
der very different conditions for forward and backward 
angles (see Sect. IIII D|l . and so this experiment was also 
subdivided into two data sets. Thus, the 16 experiments 
yield a total of 20 independent data sets for this analysis. 

The radiative corrections a ppli ed to several of the older 
experiments [^Elllilillillilllllllllilli ne- 
glected higher order terms. For the combined analysis 
of old and new experiments, the Schwinger term and the 
additional corrections for vacuum polarization contribu- 
tions from muon and quark loops have been included, 
following Eqns.(A5)-(A7) of Ref. [2J. These terms have 
very little e-dependence, and so do not have a significant 
effect on direct extractions of Ge p /Gm p from a single 
data set. However, they can modify the Q 2 -dependence 
at the 1-2% level, which has a small effect when deter- 
mining the relative normalization of the data sets. 

For some of the older experiments, there are further 
improvements that could be made to the radiative cor- 
rections, but there is not always enough information pro- 
vided to recalculate the corrections using more modern 
prescriptions. For these experiments we included only the 
terms mentioned above, which were not included in the 
earlier radiative corrections, and assume that the stated 
uncertainties for the radiative correction procedures are 
adequate to allow for the generally small differences in 
the older corrections. For a few of the earliest exper- 
iments, the quoted uncertainties for the radiative cor- 
rections were unrealistically small: total uncertainties of 
<1% or small normalization uncertainties only. To verify 
that this underestimate of the radiative correction uncer- 
tainties does not influence the final results, wc included 
a 1.5% point-to-point and a 1.5% normalization uncer- 
tainty for radiative corrections and repeated the global 
fits presented in the following sections. In most cases, 
this error was small or negligible compared to the other 
errors quoted, though for three experiments pol 1241 125| . 
the 1.5% uncertainty had a noticeable impact on cither 
the scale or point-to-point uncertainties. The additional 
uncertainties did not noticeably change the result of any 
of the fits: the extracted value of Ge p /Gm p changed by 
less than 1% for all of the fits discussed in the later sec- 
tions. Note that the results presented in this paper do 
not apply this additional uncertainty. 

Finally, we excluded the small angle data from Ref. • 
In our initial analysis, we saw a clear deviation of this 
small-angle (< 15°) data from our global fits, with or 
without the inclusion of the polarization transfer data 
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(see Sect. lIirOl or Fig.EJ. The deviat ion is due to a cor- 
rection determined to be necessary for small scattering 
angles in the analysis of NE11 that was not applied 
to the earlier Walker data [3l| . Because there is an iden- 
tified error in this data, and because it is not straight- 
forward to apply this correction to the published results, 
these small angle data (8 < 20°) are excluded from the 
analysis. 



B. Consistency Checks - Single Experiment 
Extractions 

We start by considering the values of h p Ge p /Gm p 
from published Rosenbluth extractions of the form fac- 
tors 

0, E EE E3 (Fig- El • These are the same 
experiments shown in Ref. [fjj, where the scatter is used 
to illustrate the difficulty of extracting Ge p from Rosen- 
bluth separations at high Q 2 . The data have been di- 
vided into five bins in Q 2 : 0-0.5 GeV 2 , and 1 GeV 2 bins 
above 0.5 GeV 2 . The solid line shows the weighted aver- 
age of all measurements in a given bin while the dotted 
lines show the one standard deviation range for each bin. 
As there is little Q 2 variation, averaging the values in 
small Q 2 bins should give a reasonable measure of the 
consistency of the data sets. While the average of these 
extractions is in good agreement with the global anal- 
ysis, there is a large scatter in the extracted ratios for 
Q 2 ^ 2 GeV 2 . For these measurements, \ 2 = 50.6 for 40 
degrees of freedom (45 data points minus 5 fit parame- 
ters, the mean values in each of the 5 Q 2 bins), yielding 
a x 2 per degree of freedom, \ 2 V , of 1.26 (13% confidence 
level (CL)). The agreement is worse for the higher Q 2 
data: xt — 1-63 for 17 degrees of freedom excluding data 
below Q 2 = 1.5 GeV 2 (4.9% CL). The extent of the scat- 
ter has been used to argue that the Rosenbluth extrac- 
tions do not give reproducible results. The question of 
the consistency of the data sets must be addressed be- 
fore we can draw meaningful conclusions from a global 
analysis. 

While single experiment extractions avoid uncertain- 
ties arising from the relative normalization between dif- 
ferent experiments, it should be noted that most of the 
form factor ratios shown in Fig. [21 do not correspond to 
single experiment extractions. Three of these six extrac- 
tions |3 HE E3 combine new cross section measure- 
ments with cross sections from one or more older ex- 
periments. In the extraction by Litt [24[ the new data 
are combined with results from three other experiments. 
While they give estimates of the effect of a small change 
in normalization, the quoted extractions of Ge p and Gm p 
ignore the normalization uncertainties. The extractions 
of Berger 26] and Price |27| determine normalization 
factors between their data and previous experiments by 
comparing the cross sections from the different kinemat- 
ics to the cross sections calculated assuming the dipole 
form for both Ge and Gm , in effect assuming that form 
factor scaling is valid when determining the relative nor- 
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FIG. 3: Ratio of electric to magnetic form factor from pub- 
lished Rosenbluth extractions Q, d IM IH IHIli. The 

data are binned into five Q 2 bins, and the solid lines (shaded 
regions) represent the average (one-sigma range) for the mea- 
surement in each bin. The dashed lines indicate form factor 
scaling and the fit to the polarization transfer data. 



malization factors. They do not apply any uncertainty 
associated with the determination of these normalization 
factors. 

Two of the six extractions 0, |2^| use data from sin- 
gle experiments, but use different detectors to measure 
the large and small angle scattering. Bartel [28| does 
not determine a normalization factor, but quotes a 1.5% 
relative uncertainty between the small angle and large 
angle spectrometer data. Andivahis [Tsj determined the 
relative normalization factor using data taken at identi- 
cal kinematics for the SLAC 1.6 GeV and 8 GeV spec- 
trometers. Unlike the cases where a normalization fac- 
tor between two different experiments is determined, no 
assumption about the e-dependence goes into the deter- 
mination. The uncertainty on the determination of the 
normalization factor was applied to the 1.6 GeV spec- 
trometer which provided a single, low-e point for each 
Q 2 value measured. However, the uncertainty related 
to the normalization (0.7% for Andivahis, 1.5% for Bar- 
tel) is common to all points, and will have an effect on 
Ge p /Gm p that increases approximately linearly with Q 2 
(for Q 2 <; 2 GeV 2 ). For the Andivahis measurement, 
this is roughly one-half the size of the total error, and 
so the entire data set could move up or down by roughly 
half of the total uncertainty shown in the figure. For the 
Bartel data, the high-Q 2 points all shift up or down by 
two-thirds of the total error due to a one-sigma shift in 
the normalization factor. 

Only one of the experiments Q used a single detector 
for both small and large angle scattering, and is there- 
fore free from normalization uncertainties. However, this 
is the experiment for which there was a correction for 
the small angle scattering that was not included in the 
analysis (Sect. ITlL^jl . 

To study the consistency of the Rosenbluth measure- 
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FIG. 4: Ratio of electric to magnetic form factor as extracted 
by the Rosenbluth technique, including only data sets where 
both forward and backward angle data were taken in the same 
apparatus (TablelTTl . The solid lines (shaded regions) indicate 
the average (one-sigma range) for the measurements in each 
bin. The dashed lines indicate form factor scaling and the fit 
to the polarization transfer data. The Janssens data has 19 
Q 2 points, which have been rebinned to 10 data points. 



merits without the additional uncertainty caused by com- 
bining different experiments, one must examine only 
those experiments where an adequate range of e was 
covered with a single detector in a single experiment. 
Five of the data sets from Table [i] cover an adequate 
range in e to perform a single experiment L-T separa- 
tion @> 0, H3> 0> Ell • For these experiments, we have 
used only the cross section results from the primary mea- 
surement, with the updated radiative corrections as dis- 
cussed in Sect. IIII Al 

Table|n|shows the values of G Bp /Gdi P oic, GAf p /Gdi po ic, 
and (J-pGg /Gm p from the reanalysis of the experiments 
that had adequate e coverage. Higher order terms in the 
radiative corrections have been applied to those experi- 
ments which did not include these terms, as discussed in 
Section liTTAl 

The form factors from the single experiment extrac- 
tions are given in Table ITT1 and the form factor ratios are 
plotted in Figure ^] The total \ 2 is 18-2 for 25 degrees 
of freedom (d.o.f.). If the data below Q 2 — 1.5 GeV 2 are 
excluded, X 2 = 10.3 for 9 d.o.f. (33% CL). So, while the 
published extractions of the form factors from the differ- 
ent experiments have large scatter and yield somewhat 
inconsistent results, it is in part a result of the treat- 
ment of normalization factors in these extractions. The 
raw cross sections do not show this inconsistency, and 
the true single experiment extractions are consistent, and 
agree well with the global Rosenbluth analysis. Table ITTT1 
shows the x 2 value for each of these data sets compared 
to the new global fit to the cross section data (Fig. 0, 
and the polarization transfer parameterization of Ref. [8j , 
i.e., Eqn.0] Every data set except Janssens, which is lim- 
ited to Q 2 < 1 GeV 2 ), is in significantly better agreement 



TABLE II: G_E p /Gdi po ie, GAf p /Gdi P oie, and ^pGe p /Gm p as 
extracted from the individual data sets shown in figure 2] A 
dipole mass of 0.71 GeV 2 is assumed for Gdipoic- 



Data Set 



G.Ep/Gdip GAfp/Gdip HpGe p /Gm p 



(GeV 2 



Litt 


1.499 


1 


180±0 


339 





982±0 


111 


1 


201±0 


481 


[24J 


1.998 


1 


258±0 


286 





977±0 


072 


1 


287±0 


387 




2.500 


1 


106±0 


164 


1 


010±0 


028 


1 


095±0 


192 




3.745 


1 


377±0 


256 





979±0 


038 


1 


407±0 


316 


Walker 


1.000 


1 


006±0 


072 


1 


012±0 


026 





994±0 


097 


[2] 1 


2.003 


1 


084±0 


120 


1 


027±0 


022 


1 


055±0 


139 




2.497 





944±0 


180 


1 


045±0 


022 





903±0 


191 




3.007 


1 


227±0 


145 


1 


008±0 


020 


1 


217±0 


168 


Andivahis 


1.75 





959±0 


053 


1 


049±0 


009 





913±0 


057 


[19J 2 


2.50 





863±0 


082 


1 


054±0 


008 





819±0 


084 




3.25 





868±0 


185 


1 


047±0 


015 





829±0 


188 




4.00 





890±0 


205 


1 


033±0 


015 





861±0 


210 




5.00 





578±0 


453 


1 


030±0 


016 





561±0 


448 


Berger 


0.389 





938±0 


025 





985±0 


019 





952±0 


043 


[26] 


0.584 





965±0 


019 





985±0 


009 





980±0 


027 




0.779 





950±0 


041 


1 


004±0 


012 





946±0 


051 




0.973 


1 


034±0 


058 


1 


003±0 


016 


1 


031±0 


074 




1.168 


1 


074±0 


132 


1 


022±0 


023 


1 


051±0 


152 




1.363 





907±0 


171 


1 


037±0 


022 





875±0 


182 




1.557 


1 


229±0 


265 


1 


031±0 


027 


1 


192±0 


285 




1.752 





863±0 


479 


1 


062±0 


036 





813±0 


478 


Janssens 


0.156 


1 


021±0 


028 





926±0 


027 


1 


103±0 


057 


[20] 


0.179 





962±0 


024 





959±0 


016 


1 


003±0 


039 




0.195 





973±0 


041 





999±0 


032 





974±0 


067 




0.234 


1 


020±0 


034 





939±0 


025 


1 


087±0 


061 




0.273 


1 


000±0 


039 





935±0 


019 


1 


070±0 


059 




0.292 


1 


005±0 


044 





936±0 


022 


1 


074±0 


068 




0.311 





935±0 


041 





961±0 


018 





974±0 


057 




0.389 


1 


014±0 


041 





956±0 


016 


1 


061±0 


058 




0.428 


1 


019±0 


064 





970±0 


024 


1 


051±0 


086 




0.467 





993±0 


055 





974±0 


020 


1 


020±0 


073 




0.506 


1 


023±0 


080 





954±0 


029 


1 


073±0 


113 




0.545 





984±0 


069 





983±0 


020 


1 


000±0 


087 




0.584 


1 


016±0 


103 





981±0 


030 


1 


036±0 


133 




0.623 





951±0 


085 





987±0 


020 





964±0 


103 




0.662 





869±0 


151 


1 


027±0 


031 





846±0 


169 




0.701 


1 


076±0 


100 





982±0 


021 


1 


096±0 


121 




0.740 


1 


053±0 


162 


1 


017±0 


032 


1 


036±0 


186 




0.779 





805±0 


160 


1 


035±0 


022 





778±0 


169 




0.857 





814±0 


236 


1 


083±0 


023 





751±0 


230 



1 Data below 15 degrees excluded. 

2 Using data from the 8 GeV spectrometer only. 



TABLE HI: \ 2 of single experiment extractions compared to 
Ge p /Gm from global Rosenbluth analysis and polarization 
transfer (0.6 < Q 2 < 6.0). Bartel is not included in the sum, 
as it is not taken from a single data set. 



uata oet 


data 


X (OLJ vs. 


X (OLJ vs. 




points 


global Rosenbluth 


polarization 


Litt [241 


4 


5.39 (24.9%) 


15.2 (0.43%) 


Walker [2J 


4 


5.56 (23.4%) 


20.7 (0.04%) 


Andivahis [1^ 


5 


1.44 (92%) 


13.5 (1.9%) 


Berger [2fr| 


6 


2.91 (82%) 


8.55 (20%) 


Janssens [20J 


6 


3.74 (71%) 


4.00 (68%) 


Sum 


25 


19.1 (79%) 


62.0 (0.006%) 


Bartel [28] 


8 


6.35 (61%) 


10.6 (23%) 



with the global Rosenbluth analysis. 

We could include more data sets if we also used exper- 
iments where the forward and backward angle data are 
taken with different detectors, but where the normaliza- 
tion uncertainty is taken into account. If the Bartel and 
full Andivahis data sets are included, the results are still 
consistent (x 2 =17.4 for 15 d.o.f. for Q 2 > 1.5 GeV 2 ), 
and the average form factor ratio is again slightly below 
unity for the intermediate Q 2 values. However, because 
the uncertainty related to the normalization is common 
to all Q 2 values, the Bartel data does not strongly favor 
either scaling or the polarization transfer result, and the 
reduction in the error for Andivahis when both spectrom- 
eters are combined is largely offset by the introduction of 
the correlated error. 

Because of the reduced data set and limited e range 
caused by examining only single experiment extractions, 
the uncertainties are larger than for a global analysis. 
However, this data set should be free from the additional 
systematics related to the cross-experiment normaliza- 
tion. The extracted form factor ratio is slightly below 
scaling, in good agreement with the previous global anal- 
ysis and in significant disagreement with the JLab polar- 
ization transfer results. In the next stage, we perform a 
global analysis of all of the cross section measurements 
to obtain the most precise result, and to test possible 
explanations of the discrepancy between the Rosenbluth 
and polarization measurements. 



C. Fitting Procedure and Results 

For the global fits to the cross section data, the form 
factors arc parameterized using the same form as Refs. |9j, 



G E {q),G M (q)/n P = 1 / [1 + Piq + P2q 2 - 



■PNq 



AT] 



(6) 



where q = y^, and TV varies between 4 and 8; N=5 is 
adequate for a good fit. Because we wish to focus on the 
discrepancy at intermediate Q 2 values, we exclude data 



below Q 2 = 0.6 GeV 2 , where the two techniques are in 
good agreement, and above Q 2 = 8 GeV 2 , where the data 
do not allow for a Rosenbluth separation. Several param- 
eterizations were tried, and this form was chosen because 
it provides enough flexibility to fit the Q 2 -dependence of 
the form factors. This parameterization has odd powers 
of q, and so will not have the correct Q 2 — > behav- 
ior, but this is not relevant for this analysis, as we are 
focusing on the intermediate Q 2 range. 




2 3 4 
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FIG. 5: Ratio of electric to magnetic form factor from the 
new global fit (solid line), the Bosted fit (l^t (dashed line), 
and the fit to the polarization data (dotted line) Q]. 

In addition to the 2N parameters for the form factors, 
we also fit a normalization factor for each of the data 
sets. After partitioning the data sets for experiments 
that make multiple independent measurements, there are 
20 normalization constants for the 16 experiments, as 
described in Sect. IIII Al The fit parameters are allowed 
to vary, in order to minimize the x 2 , the total x 2 for the 
fit to the cross section data: 



xl 



£ 

i=l 



+ 



N cxp 



(v 3 i) 2 



(7) 



where <7i and dai are the cross section and error for each 
of the N„ data points, r\j is the fitted normalization fac- 
tor for the j'th data set, and drjj is the normalization 
uncertainty for that data set. From the fit, we obtain 
the values of the fit parameters, pi (Eqn. EJ) , for both the 
electric and magnetic form factors, as well as the nor- 
malization factors for each of the 20 data sets. Figure El 
shows the result of our global fit along with the fit of 
Bosted 0| to the previous global analysis and the fit to 
the polarization transfer data. While the modifications 
described in Sect. IIII Al do decrease the h p Ge p /Gm p for 
Q 2 > 2 GeV 2 , the new fit still gives a result that is well 
above the polarization transfer result at these momen- 
tum transfers. The x 2 value for this fit , obtained from 
Eqn.[7| is 162.0 for 198 d.o.f. (i.e., xl = 0.818). 
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FIG. 6: Ratio of electric to magnetic form factor as deter- 
mined from direct L-T separations, using the normalizations 
determined from the global fit. The dotted line is the ratio 
from the global fit to Ge p and Gm p , while the dashed line is 
the fit to the polarization measurements. 



In order to estimate the uncertainty of the fit as a func- 
tion of Q 2 , we perform direct L-T separations wherever 
there are enough data points in a small range of Q 2 . The 
normalization factors from the global fit are used to scale 
each data set, and the Q 2 -dependence of the fit is used to 
scale each point to the central Q 2 value. Figure [S] shows 
the extracted ratios and uncertainties from these direct 
L-T separations at 26 Q 2 points up to 6 GeV 2 (15 above 
Q 2 =0.6 GeV 2 ). These Q 2 values are selected by requir- 
ing that five or more e values are in each Q 2 bin, that the 
e-range of the data is at least 0.6, and that the correction 
required to scale the data points to the central Q 2 value 
never exceeds 2%. Due to this constraint on the correc- 
tion, the resulting form factor ratio is independent of the 
model used to scale the cross sections for any reasonable 
model of the form factors. 

To bring the data into agreement with the polarization 
transfer measurements, there would have to be a signif- 
icant e-dependent error in the cross sections. Assuming 
that all of the data sets have such an error, and that it 
is linear in e, it would have to introduce an e-dependence 
of 5-6%, nearly independent of Q 2 , for Q 2 <; 1.0 GeV 2 . 
This error would have to be even larger if it affected only 
some of the data sets, or if all of the uncertainty came at 
very large or very small values of e, e.g., if the error only 
occurred at very small angles or very low energies. In ad- 
dition, a repeat of the above fit, excluding data at high 
(e > 0.8) or low (e < 0.3) values of e did not have sig- 
nificant impact on the overall fit, changing the extracted 
values of G Ep by <5% for Q 2 values below 4 GeV 2 . At 
higher Q 2 values, the cuts on e led to larger changes in 
Ge , but the reduction in the e range did not allow for 
a precise extraction of Ge , and the changes were small 



compared to the final uncertainty. 



D. Consistency Checks 

As a first check of the consistency of the data sets, we 
examine the contribution to the total \ 2 coming from 
each data set. No data set had an excessively large x 2 
value, only five of the data sets have x 2 > 1 and all five 
of these have reasonable confidence levels (none below 
10%). In addition, no individual cross section value had 
excessively large (> 3er) deviations from the fit. However, 
the fact that the fit gives such a low xt value, combined 
with the fact that most of the individual data sets have 
X 2 < 1, indicates that several data sets may have overly 
conservative estimates for the uncertainties. Therefore, 
the total x 2 value of the global fits cannot be viewed as 
an absolute measure of goodness of fit, and we will focus 
on the change in x 2 between different fits when using the 
same data sets as a measure of the relative goodness of 
fit. 

While these statistical measures help us locate individ- 
ual data sets that are inconsistent with the bulk of the 
data, they are not always enough to detect systematic 
errors in the data, especially when the error estimates 
are somewhat conservative. It is quite possible that these 
statistical measures will overlook systematic errors which 
are small compared to the individual uncertainties, but 
still large enough to modify the small e-dependence ex- 
tracted. Thus, we would also like to look for any sys- 
tematic trends in the data sets which might modify the 
result. 

For each data set we compare the measured cross sec- 
tions to the fit and look for systematic deviations from 
the global fit as a function of Q 2 , 8, E' e or e. A sys- 
tematic e-dependence of just 1-2% may not show up in 
the total x 2 °f t ne data set (especially if the errors are 
overly conservative or are noticeably larger than 1-2%), 
but could lead to a noticeable change in the electric form 
factor at high Q 2 values. In the preliminary version of 
the global fit, it was observed that the Walker data @ 
showed a clear deviation from the global fit at small an- 
gles (Fig. 0|, even though the x 2 value for this data set 
was quite reasonable. The data below 20° was then re- 
moved from later fits. 

Similarly, the Goitein data [25| showed a systematic 
deviation from the initial global fit. The ratio of cross 
section from this data set to that obtained from the global 
fit was systematically higher by <; 5% for the higher Q 2 
values. The original publication listed several differences 
in the larger Q 2 data: different collimation, additional 
cuts, and several corrections which were negligible for the 
low Q 2 data were quite large for these points. Because 
of the difference in running conditions between low and 
high Q 2 , it was decided to break up this data set into 
two subsets, each with its own normalization factor, and 
to increase the normalization uncertainty from 2.8% to 
3.8% for the large Q 2 data due to the larger and less well 
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FIG. 7: Ratio of cross section to global fit for Walker data |2|, 
showing a systematic deviation at low scattering angles. 

understood corrections for this data set. The high Q 2 
data could also have been excluded altogether, but while 
the conditions for this subset of the data were different, 
there was no clear evidence of any specific problem. 

As another check for systematic trends in the data, we 
examine the direct L-T separations performed at several 
values of Q 2 . For the separations shown in Fig. El we plot 
the cross sections versus e to look for systematic devia- 
tions from the expected linearity and to look for data sets 
which have systematic differences in their e-dependence. 
After separating the Goitein data into a two subsets 
the only data set which stands out is the data from 
Janssens j2ij, which has some deviations from a linear 
e-dependence. Sometimes the higher e points are above 
the extracted (linear) e-dependence, sometimes they are 
below, and sometimes they just show non-statistical scat- 
ter. Because there are no systematic trends to these devi- 
ations, and because there was no indication of problems 
in the original publication, we do not exclude this data 
set from the analysis. The effects of excluding specific 
data sets was tested separately, and the results are pre- 
sented in the following section. 



E. Bad Data Sets? 

Because most of the high Q 2 Rosenbluth data comes 
from a limited set of experiments, it is possible that a sin- 
gle data set may have an error that strongly biases the 
global analysis. The global fit described in Sect. lIII HI was 
repeated 20 times, with a different data set excluded each 
time. The exclusion of a single data set generally had lit- 
tle effect on the global fit, although there were a few data 
sets whose exclusion lead to a noticeable increase (up to 
~ 10% for Q 2 > 2 GeV 2 ) or a noticeable decrease in 
H P G Ep /G Mp (from 5% to 15% for 3 < Q 2 < 5 GeV 2 ). 
Even excluding two or three data sets together gener- 
ally had little effect on the result. The three data sets 
whose exclusion lead to the greatest decrease in the ratio 
were identified, and the fit was repeated with all three 
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FIG. 8: Ratio of electric to magnetic form factor from the 
global fit of all of the data sets (solid curve), and the lowest 
ratio found by excluding up to three of the data sets (dashed 
curve). Also shown is the fit to the polarization data (dotted 
line) . 

excluded. Figure |H1 shows the result of the full global fit, 
and the global fit with all three of these data sets re- 
moved. There is little change for Q 2 < 3 GeV 2 , where 
both techniques have high precision, while there is a sig- 
nificant decrease for Q 2 > 3 GeV 2 , where the ratio is 
constrained by a very limited set of experiments and is 
very poorly constrained when the three data sets are re- 
moved. Also, while the ratio at large Q 2 is significantly 
decreased, it is still well above the polarization transfer 
result. 



F. Consistency of Global Fit and Polarization 
Transfer Results 

The new global fit is clearly in disagreement with the 
polarization transfer results, and shows no indications of 
inconsistency between the data sets or bias due to in- 
clusion of erroneous results. One remaining possibility 
is that it is the fitting procedure itself, rather than the 
data, that leads to this discrepancy. While we include 
the relative normalization factors for the experiments in 
the analysis and obtain best-fit values for the normaliza- 
tions, it is possible that a small change in the normal- 
ization factors for some or all of the experiments would 
significantly improve the agreement with the polarization 
transfer data without dramatically increasing the overall 
X 2 of the fit. 

To test this possibility, we performed a constrained fit, 
using the same data as in the previous fit, but fitting only 
Gm and the normalization factors, while fixing the ratio 
H p Ge /Gm to the fit of the polarization transfer data 
(Eqn. . The extracted magnetic form factor is roughly 
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2% higher over the entire Q 2 range (Q 2 J> 1 GeV 2 ), due 
to the reduction in strength from the electric contribu- 
tion. The constrained fit has 5 more degrees of freedom 
(the parameters that were used to fit Ge p ) while the to- 
tal x 2 for the fit increases by 40.2 (from x 2 =162.0 for 198 
d.o.f. to x 2 =202.2 for 203 d.o.f.). 

Even though the normalization factors are adjusted 
in order to best reproduce the polarization transfer re- 
sults, the direct L-T separation using these normaliza- 
tion factors is systematically high for all Q 2 values, as 
seen in figure OH So not only is this fit significantly 
worse, the normalization factors derived from constrain- 
ing the Ge p /Gm p ratio to match polarization transfer 
data still do not fully explain the magnitude of the falloff 
of Ge p /Gm p - The effect is only slightly less using the 
combined fit to all four recoil polarization data sets. With 
the ratio forced to n P G Ep /G Mp = 1 - 0.135(Q 2 - 0.24) 
(Eqn. EJ), the total x 2 is 197.3, an increase of 35.3 in the 
total \ 2 , with 5 extra degrees of freedom. 

Because of the conservative error estimates on some of 
the data sets, the reduced x 2 is still less than one, even af- 
ter this noticeable increase. Therefore, the absolute value 
of xt cannot be directly used to measure goodness of fit. 
To get an idea of the relative goodness of fit, one can 
scale the overall uncertainties on the cross sections such 
that xt f° r th e unconstrained fit is approximately one. 
This means reducing the cross section uncertainties by a 
factor of 0.905, leading to x 2= l f° r the unconstrained fit 
to the cross section data, x 2 = l-22 (1.9% CL) for the fit 
constrained to Eqn.H and x 2 =1.19 (3.5% CL) for the fit 
constrained to Eqn.0 Assuming that the unconstrained 
fit should yield x 2 = 1 is arbitrary, but it is a reasonable 
staring point since we expect x 2 ~ 1 if the estimated 
uncertainties are correct and if our fitting function can 
accurately reproduce the data. One would expect x 2 to 
be even higher, and the confidence levels for the con- 
strained fits to be even lower, if the fitting function does 
not adequately reproduce the data or if there are any 
inconsistencies in the cross sections. So it is likely that 
these confidence levels are upper limits of the consistency 
between the cross section data and the parameterizations 
of the polarization transfer data. 

Forcing the fit to match the parameterization of the 
polarization transfer data gives too much weight to polar- 
ization data, as it neglects the uncertainties in the polar- 
ization measurements. To avoid this, we also performed 
a combined fit, treating the cross section and polarization 
transfer measurements on an equal footing. We repeated 
the procedure described in Section IIII CI but with the 
inclusion of the polarization transfer ratios as additional 
data points, and with a systematic offset included for 
each data set (H HI > as in the fit from section III Bl 
(Eqn. [SJ . The x 2 f° r this combined fit is the contribu- 
tion from the cross section measurements (Eqn. 0) plus 
the additional contribution from the polarization transfer 
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FIG. 9: Ratio of electric to magnetic form factor as deter- 
mined from direct L-T separations, using the normalization 
factors determined by forcing the ratio to reproduce the polar- 
ization transfer fit (dashed line). Although the normalization 
factors were fitted to reproduce the polarization transfer re- 
sults, the direct L-F separation is systematically higher than 
the fit at large Q 2 values. 



ratio measurements: 

n r 



X 



/Co 



E 

k=l 



[Rk — Rat) 2 



{dR s 



E 

i=i 



(AO 2 

{dRsys) 2 



(8) 



where R = Ijl p Ge p /Gm p , rf-Rstat and dR sys are the statis- 
tical and systematic uncertainties in R, Nr is the total 
number of polarization transfer measurements of R, A; 
is the offset for each data set, and N exp is the number of 
polarization transfer data sets. Because the polarization 
transfer data are included in the fit, the normalization 
factors for the cross section measurements are adjusted 
to give consistency between different cross section data 
sets, as well as consistency with the polarization mea- 
surements of fipGE /Gm ! much as they are in the con- 
strained fit. 

The ratio of Ge /Gm for this combined fit is sys- 
tematically higher then the polarization transfer results 
(Fig.llUp. As with the constrained fit, the direct L-T sep- 
arations using the normalization factors from this global 
fit are systematically higher than the actual ratio ob- 
tained in the fit, and the quality of fit is significantly 
reduced when the polarization transfer data is included: 
X 2 = 215.8 for 218 d.o.f., an increase in x 2 of 53.8 for the 
additional 20 data points. 

While the combined fit to cross section and polariza- 
tion transfer data yields a ratio that is close to the polar- 
ization measurements, this does not imply that the data 
sets are yielding consistent results, just that the polariza- 
tion transfer results dominate the fit. The goodness of fit 
is significantly worse when the polarization transfer data 
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FIG. 10: Ratio of electric to magnetic form factor as de- 
termined from combined fit of cross section and polarization 
transfer data. The dotted line is the result of the fit, and the 
circles are the results from the direct L-T separations using 
the normalization factors determined from the global fit. 



is added to the cross section measurements. Because of 
the conservative errors, it is difficult to estimate the rel- 
ative goodness-of-fit for the combined fit. We can scale 
the uncertainties by a factor of 0.905 to force xt — 1 
for the fit to cross sections, yielding xt = 1-21 (1.9% 
CL) for the combined fit to the cross sections and po- 
larization transfer. However, as discussed for the con- 
strained fit, this is likely to be an upper limit on the 
consistency of the two fits. The comparison of single ex- 
periments L-T separations to the polarization transfer 
fit, HpG Ep /G Mp = 1 - 0.13(Q 2 - 0.04), yields a confi- 
dence level of 0.006% (Table EQJ. To take into account 
the uncertainties of the polarization transfer measure- 
ments, we compare the single experiment L-T separa- 
tions to the fit with the slope parameter decreased by 
one sigma, i.e., fi P G Ep /G Mp = 1 - 0.125(Q 2 - 0.04) [jf. 
This gives x 2 = 57.8 for 25 degrees of freedom, or a 0.02% 
confidence level. 

We can also use a recently-proposed test to determine 
the consistency of the data from the two techniques, fol- 
lowing the prescription of Refs. [13, HH- They define 
a test called the "parameter goodness-of-fit" (PG) mea- 
sure, which is designed test the consistency of indepen- 
dent data sets sharing a common set of parameters. This 
approach is significantly less sensitive to over- or under- 
estimates of the uncertainties of the individual data sets. 
The corresponding x 2 for this test, x 2 m in; * s the increase 
in x 2 f° r a combined fit, relative to fits of the individ- 
ual data sets. The effective number of degrees of free- 
dom is the number of parameters in common between 
the two data sets, which in this case are the parame- 
ters that go into the ratio. So, for comparison of the 
cross section and polarization transfer results, X 2 min i s 



the difference between Xmin f° r the combined fit and 
the values of Xmin f° r the independent fits of the cross 
section data and of the polarization transfer data, i.e., 
X 2 min = 215.8 - 162.0 - 26.5 = 27.3. Note that the 
X 2 for the fit to polarization transfer data alone is not 
the result presented in section 111 Bl because the com- 
bined fit includes only the polarization transfer data from 
JLab |6J, [?], |8| . The two types of data both constrain the 
ratio, effectively 5 parameters, and so the PG goodness- 
of-fit measure corresponds to x 2 min = 27.3 for 5 degrees 
of freedom, or a 0.005% CL. This measure has two advan- 
tages over the confidence levels calculated based on the 
direct comparison of the single experiment extractions to 
the polarization transfer fits. It includes the entire cross 
section data set, rather than just the five single exper- 
iments extractions in Table II I II and it accounts for the 
uncertainties of the both the cross section and polariza- 
tion transfer measurements. 

In the end, both approaches yield a confidence level of 
•C 1%, showing that the data sets are clearly inconsis- 
tent. If it is an error in the cross section data, it would 
take an e-dependent correction of at least 5-6% to make 
the results from the two techniques consistent. Until the 
source of the inconsistency between the two data sets is 
determined, we do not know how to combine the mea- 
surements in order to obtain a reliable extraction of the 
form factors. 



IV. DISCUSSION 

We have presented a new extraction of the proton elec- 
tromagnetic form factors based on a global analysis of 
elastic cross section data at moderate to high Q 2 . The 
extracted value for Ge p /Gm is slightly lower than in 
the previous global analysis 0, 0] , but is still well above 
the values determined by polarization transfer measure- 
ments. We have demonstrated that the discrepancy be- 
tween the global Rosenbluth analysis and the polariza- 
tion transfer results is not merely the result of the in- 
clusion of one or two bad data sets. While modifying 
the normalization factors of the different experiments can 
lead to a significant reduction in the extracted Ge p /Gm p , 
they yield significantly lower quality fits, and still do not 
fully resolve the discrepancy between the two techniques. 
Finally, we have demonstrated that the apparent incon- 
sistency between single experiment extractions is not due 
to any inconsistency in the data themselves, but due to 
improper treatment of normalization uncertainties when 
combining data from multiple measurements. The values 
of Ge p /Gm p from single experiment extractions, which 
avoid the normalization uncertainties involved in a com- 
bined analysis, are self-consistent and in good agreement 
with the result of the global analysis. 

This indicates that there is a more fundamental reason 
for the discrepancy, such as an intrinsic problem with ei- 
ther the Rosenbluth or polarization transfer techniques, 
or an error in the cross section or polarization transfer 
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measurements. If the error is in the cross section mea- 
surements, it must be a systematic problem that yields 
a similar e-dependence in a large set of these measure- 
ments: a 5-6% linear e-dependence in all cross section 
measurements above Q 2 ss 1 GeV 2 , larger if only some 
of the data are affected, or if the error has significant de- 
viations from a linear e-dependence. Since the the cross 
sections are necessary to extract the absolute value of 
Ge and Gm even with very precise measurements of 
Ge I Gm , an unknown systematic error in the cross sec- 
tion measurements implies unknown systematic errors in 
the values of the form factors. 

Until this discrepancy is understood, it is premature 
to dismiss the Rosenbluth extractions of the form fac- 
tors. There is a significant difference in the values of Ge p 
extracted by these two techniques, and smaller (~3%) 
differences in the extracted values of Gm p - In addition 
to the impact this uncertainty has on the state of our 
knowledge of the proton structure, it also affects other 
measurements which rely on the proton electromagnetic 
form factors as input, or measurements where there are 
significant corrections due to the radiative tail of the elas- 
tic peak. 

While the cross sections extracted using the new polar- 
ization transfer ratio measurements are typically within 
a few percent of previous parameterizations, the large 
change in the ratio of the form factors means that 
the changes in extracted cross sections are strongly e- 
dependent. Thus, it will have a significant impact on 
experiments that try to examine a small e-dependence. 
For example, there will be a large effect on Rosenbluth 
separation measurements for A(e, e') (e.g., Coulomb sum 
rule measurements at large q) or A(e, e'p) [e.g., > 50% 
difference for recent results for separated structure func- 
tions in nuclei [SHH)- 

Experiments which need only the unseparated elastic 
cross sections are less sensitive to these uncertainties. 
Two such examples are the extraction of the axial form 
factors from neutrino scattering measurements, and ex- 
tractions of nucleon spectral functions in nuclei, which 
use the proton electromagnetic form factors (or elastic 
cross sections) as input. An analysis of the neutrino mea- 
surements 36] shows that the extracted axial form factor 
is essentially identical using the Rosenbluth or polariza- 
tion transfer parameterizations of the form factors. The 
difference also has relatively little effect on the extraction 
of the spectral function in unseparated A(e, e'p) measure- 
ments. This is not surprising, as the polarization mea- 
surements only determine the relative contributions from 
Ge and Gm , and the overall size of the form factors 
is still determined by fitting the cross section measure- 
ments. With the constraint from the polarization mea- 
surements included, the combined fit can reproduce the 
measured cross sections at the few percent level. How- 
ever, while fits with and without polarization transfer 
data can reproduce the measured cross sections at the 
few percent level, the discrepancy in the ratio of form 
factors implies that some large set of these cross sections 



may be wrong by 5% or more. Until the source of the dis- 
crepancy is identified, there is no way to be sure which of 
these cross sections are incorrect, or how well these cross 
sections are measured. 

Finally, it is important to note that the while a consis- 
tent extraction of Ge and Gm yields only small differ- 
ences in the elastic cross sections, this is not always the 
case when combining extractions of Ge and Gm from 
different analyses. Using the polarization transfer data 
to constrain Ge /Gm and cross section measurements 
to determine the size of Ge p and Gm p , as presented here 
or in another recent analysis , yields cross sections that 
are nearly identical to those from the cross section anal- 
ysis only, except for a change of approximately 5% in 
their e-dependence. However, using Gm as extracted 
from Rosenbluth extractions, combined with Ge /Gm 
from polarization transfer measurements can yield cross 
sections which are significantly lower than the measured 
cross sections. For example, if one takes the magnetic 
form factor parameterization from Bosted [l5j , but calcu- 
lates the electric form factor using the polarization trans- 
fer ratios, the resulting cross sections are lower by 4-10% 
over a very large Q 2 range (0.1-15 GeV 2 ), compared to 
using both form factors from the Bosted parameteriza- 
tion. Thus, it is extremely important to use a consis- 
tently determined set of form factors when examining 
the difference between the Rosenbluth and polarization 
transfer data, or when parameterizing the clastic cross 
section for comparison to other data or calculations. 

There have been two more recent L-T separations from 
Hall C at JLab, which were not included in this analy- 
sis. One took points at Q 2 =0.64 and 1.81 GeV 2 |35J . 
both of which were in agreement with the new global 
Rosenbluth analysis presented here (Fig. A more 
extensive set of Rosenbluth measurements was taken as 
part of JLab experiment E94-110 [33]. The experiment 
measured R — o~l /&t m the resonance region, but also 
took elastic data. This allowed them to perform several 
Rosenbluth separations for 0.5 < Q 2 < 5.5 GeV 2 [38| . 
These results are also in excellent agreement with the new 
global fit presented here. While these newer JLab mea- 
surements have not been included in the new fit, the fact 
the single experiment extractions from these measure- 
ments are in good agreement with the global fit indicates 
that they would not significantly change the final fit, al- 
though their inclusion would decrease the uncertainties 
in the global analysis, and thus increase the significance 
of the discrepancy with the polarization transfer results. 



V. SUMMARY AND FUTURE OUTLOOK 

A careful analysis of Rosenbluth extractions has been 
performed to test the consistency of the worlds body of 
elastic cross section measurements, and to test explana- 
tions for the discrepancy between polarization transfer 
and Rosenbluth extractions of the proton form factors. 
We find no inconsistency in the cross section data sets, 
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and cannot remove the discrepancy via modifications of 
the relative normalization of different data sets or the 
exclusion of individual measurements. 

This discrepancy indicates a fundamental problem in 
one of the two techniques, or a significant error in either 
the polarization transfer or cross section measurements. 
An error in the polarization data would imply a large er- 
ror in the electric form factor extracted from a combined 
analysis, and may have consequences for other recoil po- 
larization measurements. An error in the cross section 
data would have to introduce an e-dependence of > 5% 
for Q 2 > 1 GeV 2 , implying an error in both the elec- 
tric and magnetic form factors However, while the un- 
certainty in the form factors is smaller in this case, the 
error in the cross section measurements will also lead to 
uncertainties in other measurements which require the 
elastic cross section as input. Thus, even if it is demon- 
strated that the polarization transfer measurements arc 
correct, it is necessary to determine the source of the 
discrepancy in order to have confidence in our knowledge 
of the elastic cross sections, and any other measurement 
which relies on this knowledge. 

Future results from JLab should significantly improve 
the situation. JLab experiment E01-001 ran in 2002, and 
performed a high precision Rosenbluth separation using 
a modified experimental technique |39 . |40| . This experi- 
ment should be able to clearly differentiate between the 



ratio of Ge p /Gm p as seen in previous Rosenbluth mea- 
surements and in the polarization transfer results, while 
being significantly less sensitive to the types of system- 
atics that are the dominant sources of uncertainties in 
previous results. In addition, there is an approved ex- 
periment to extend the polarization transfer measure- 
ments to higher Q 2 |4l| using a different spectrometer 
from the previous measurements (Q 2 > 1 GeV 2 ), which 
all used the High Resolution Spectrometer in Hall A at 
JLab. This will provide the first independent check of the 
Hall A large Q 2 polarization transfer experiments, whose 
systematics are dominated by the spin-precession in the 
spectrometer. These tests will provide crucial informa- 
tion to help explain the discrepancy in the measurements 
of the proton form factors. The resolution of the discrep- 
ancy will significantly improve the state of knowledge of 
the proton form factors, as well as determining if other 
other measurements utilizing the Rosenbluth or polariza- 
tion transfer techniques are affected. 
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